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ABSTRACT 


The analysis of non-linear systems has frequently been a major 
problem to the engineer. The solution of system equations often 
requires either a computer or relatively complex numerical tech- 
niques. A straight-forward, relatively simple method is proposed 
herein which permits the engineer to satisfactorily approximate the 


exact solution to certain non-linear differential equations. 
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I. INTRODUCTION 


Systems described by non-linear differential equations often 
must be analyzed on a computer, using some numerical algorithm if 
the computer is of the digital type. The system variables are 
evaluated at specific points, and perhaps graphs drawn of these 
points. To discover the actual time solution to the system equations 
it is usually necessary to fit these points to some other numerical 
algorithm which will provide a curve fitting polynomial form of the 
exact solution. Such methods not only are not completely accurate 
but also are restrictive in that they require computer access; also, 
depending upon the complexity of the particular system, computer 
solution may require excessive time. These restrictions can sometimes 
be avoided if there is available a mathematical method to analyze the 
non-linear system. 

Often, engineers are not as interested in exact solutions as 
they are in answering the questions: 

1. What are the initial and steady-state values of the system 
variables? 

2. What are the transient characteristics of the system? If the 
engineer can answer these questions he has, in effect, analyzed the 
system. The general nature of these questions suggests that the exact 
solution is not necessary if an approximate solution can be found which 


has essentially the same transient and steady-state responses. 


TI. DEVELOPMENT OF THE APPROXIMATION METHOD 


Since the objective of an approximation method is to make the 
solution to non-linear equations readily obtainable, a starting 
point in finding such a method is to try to manipulate the equations 
into a linear, or at least less non-linear form. A method which has 
characteristics that can accomplish this purpose is the Laplace 


Transform equivalent to time equations. 


A. DERIVATION OF THE SERIES FORM OF THE LAPLACE TRANSFORM 

Baycura [ref. 1] proposed that the transform of certain non-linear 
terms could be found by using the series form of the Laplace Transform. 

Recall that X(s) = [ x(eye"Stat (II-1) 
Expanding the oigive liek) 2400 by integrating by parts, 

fered 

[xerStat = xX [estat - x { [ e°Statae + x [ff erstatatat-... 

aiid hein x, ‘its derivatives, and the “vena s at the extreme 


points, it can be shown that the expansion becomes 


t{x} = x(0) 4 x(0) , x) x(0) ex (II-2) 
S S 
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Using this series form of the Laplace Transform, Brady [ref. 2] found 


the transform of x"(t). 


PRs as (11-3) 
Using similar methods a calculable final value term can be used as one 
of the extremes in the evaluation of the Laplace Integral to obtain a 
linearized series approximation to the Laplace Transform of non-linear 


terms. 








B. OBTAINING THE FINAL VALUE 

The method used to derive the series form of the Laplace Transform 
was to expand the Laplace Integral by integrating by parts and then 
evaluate each term at the extreme points, a necessary condition being 


that xe>° 


+0Oas t —+o. The upper limit represents the steady-state, 
and for certain non-linear differential equations the steady-state, or 
final value, can be found. 

For example, consider the non-linear differential equation 


yo ae xe = 1, x(0) = 0 (11-4) 


In the steady-state, x = 0, and 


2 
O + Xe 


Solving this quadratic equation produces the roots: 


7 Xe = 1, where Xe is the final value of x. 


Xe = G15, -18605 


C. APPROXIMATING THE LAPLACE TRANSFORM OF THE NON-LINEAR TERM 
After integrating by parts, the expanded form of the Laplace Trans- 
form of x2 is 


fix?) = x [xe Stat - x [ [restarat + 8/ I fre" Fataede--- (II-5) 
oO oO oo 


° 


Each term contains the term xe Stat which is, by definition, the 

oO 
Laplace Transform, X(s). To obtain the series approximation each term 
in equation (II-5) is evaluated at its extreme points. In this manner 


the first term becomes 


x [xe at = (X - Xq) *X(s) 


Since integrating in the time domain corresponds to multiplying by L 


in the s domain, the second term becomes 


- X,)X 
i [TaeStaeds =~ Gig = ig) nfapae = - StS 


oO 


It follows that similar evaluation of al] the terms in the expansion 
will result in an infinite series of the form: 
(Xe - Xp) K(s) (Bp - Ky) X(s) 


(X¢ - Xq)X(s) 7 — = s , 


Six?) 


(Xe - X) (Xe = Xp) 
= [ (Xp - %) --t—2.+-t,+ > rae (TT=6 


S 

D. THE METHOD OF APPROXIMATION USING SUCCESSIVE TERMS IN THE INFINITE 

SERLES 

The terms in the infinite series form given in equation (II-6) can 
be evaluated by taking derivatives of the original differential equation. 
Including more and more terms in the series will result in successive 
approximations to the transform of the non-linear term. Evaluation of 
these series terms is illustrated using equation (II-1). With the 
given conditions, x, = 0 and Xe = 0, 
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X+xX+tx =] Xe = .615, Xo = ] 
Xx +x + 2xx = 0 Xo = -1, Xe = 0 
+H + xX + 2k% = 0 KG = -1, KE = 0, etes. 


E. USE OF THE FINAL VALUE IN THE APPROXIMATE TRANSFORM 

Using the series form of the transform of x2, a final value term 
can be inserted to yield a linear equation in X(s) that is an approxi- 
mation to the transform of 2. 

For example, consider a first approximation using only the first 


term of the series. 


Six2} = X(s)(X¢ - X) = XeX(s) 








Transforming equation (11-4) term by term, 
tx?) = §X(s) - Xq + X(s) + s6NI5X(S) =< 


ieiieet 1615) = . 


we 
MENS == + 1.615) 


This equation can be expanded in partial fractions and the inverse of 


the s-domain solution found. 


i> 


; _A B 
SS) = eee S| «Ss CERIO TS 


2 ee ee 2 ml __ 
| memos 's-0 Oem = |'s = -1.c15uaue-°!® 


Pain ae 
NS) Ss aadiars, 615 (11-7) 


The time solution corresponding to equation (II-7) is 


x(t) = .618 - .618 e /O15t 

This approximate solution has an initial value of zero and rises 
exponentially to a final value of .618. These values compare quite 
favorably with the given initial condition, XQ = 0, and calculated 
final value, Xe = .615. A comparison between the approximate solution 
and the computer solution found on the SDS 9300 digital computer is 
Shown in Table II-1] and Figure II-1. The comparison shows a close 
agreement between the two transient responses. 

A second approximation to the transform would include the second 


term in the series form. 
2. \( (X¢ > ‘0 | 
ie = Xs) | (X¢ - XQ) secre maa 


Substituting the values calculated for the terms in the series the 


2 


second approximation to the transform of x~ becomes SHS} = X(s)(.615-z). 


The effect of including derivative terms in the approximation is 


discussed in later sections. 
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. OOZAIG 
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44822 
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-49509 
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sO le 
.6178] 
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TABLE II-] 


ERROR 


. 00000 
.00270 
.00853 
.02038 
O22 
.02650 
.03152 
.03421 
.03560 
.03586 
0882 
.01479 
.00387 
.00075 
.0002] 
.00000 
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a2 Bo 4 


time(sec. ) 


FIGURE II-] 
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III. APPLICATION OF THE FINAL VALUE METHOD OF APPROXIMATION 


The use of the final value method of approximating the solution 
to non-linear system equations can now be demonstrated on two equations 


containing first order non-linear terms. 


A. APPLICATION TO A NON-HOMOGENEOUS NON-LINEAR EQUATION 
Consider the case of vertical fall with air resistance. The 


governing equation of the system is 
MY = Mg - kv*, v(0) = 0 (III-1) 


where M is the mass of the falling object, g is the acceleration due 
to gravity, and k is a constant relating to the retarding force, which 
varies as the square of the velocity. Evaluating the system equation 


in the steady-state where v = 0, the final value of the velocity is 


determined. 
a 
Oigsts KVe Mg 
ee 
Ve = Mg/k 
Vee \Mg/k 


The series approximation to the Laplace Transform of ve as given by 
equation (II-6) is 

SANPT = V(s) | (ve - Vo) = pi ae] 
Having obtained the final value and knowing the series form of the 
approximate transform of Vv", Successive approximations can be made to 
the solution of the system equation, adding successive terms with each 


approximation. 
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1. First Approximation. 


Using the first term in the series form only, 


£iv*} = veVis) = \Wg/k V(s) 
Transforming the system equation term by term, 
LIM = Mg - kv°} = sV(s) + (k/M)veV(s) - g/s = 0 
V(s) = sis + ve/M 


Expanding in partial fractions, 


Vs s(s + kv /M) - 


Mg/kv - Mg/kv ¢ 
Vis) = —S—— = Fey 


Taking the inverse transform, and substituting the equation vee \Mg/k, 


B 
+ 
Ss + kv -/M 


wn [> 


the time solution is found. 


a v(t-e° (YAM (11-2) 


The exact solution can be obtained by separation of variables: 
¥ 2 
dv/dt = g - kv-/M 


Ws dt 
g - kv -/M 


kM eee 


Mg/k - v 


t = kM(1/2 Mg/k)In wae 


Mi sartesah/ 


2 
e° gk/M = 
Ve V 


-2 
— 1, ae mt -3) 
(1 +e“ {gk/M't) 
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As in the exact solution the velocity given in the approximate solution 
starts at zero and increases exponentially to Ve in the steady-state. 
The basic difference between the two solutions is the speed with which 
the velocity reaches the steady-state value. A comparison between the 
two solutions is shown in Table III-] and Figure III-1], with Ve = 36.4 
and \gk/M = .27. 
2. Second Approximation. 
In this Second Approximation the first two terms in the series 


form are used and the approximate Laplace Transform becomes: 


(Vv. - V,) 
ive} + v(s) |v, - aia (III-4) 


The added terms are evaluated remembering that the necessary condition 
for the final value method of approximation is that Ve = 0. Vo can be 


calculated from the system equation. 


Vo > gM/M - Ok 


Substituting this value into equation (III-4), 


Ziv} = V(s) [V- + g/s] 


After transforming the system equation term by term, a partial fraction 


expansion is used to find the s domain characteristic equation. 


sV(s) + V(s)(k/M)(v_ + g/s) = g/s 


(Ss + ky 
V(s) 


(8) = SERRA SF GEM 


Expanding in partial fractions and taking the inverse transform, the time 


f 4 kal) 5 o/s 


solution is found. With L = k/M, 


ee 
(s + \Lg/4 - j \BLg/4)(s + \Lg/4 + j \Blg/4) 


_. ee A eee Ax 
(s + {Lg/4 - j Y3Lg/4) (s + WLg/4 + j YBLg/4) 


a! 
jBlg Bbq \3Lg 


v(t) = \W73 \ezL ec” NE9/4 tein Boze t (III-5) 


It is apparent that v(t) has a steady-state value of zero. Also, this 
approximation has the transient characteristic of a damped sinusoidal, 
while the exact solution was essentially exponential with no oscil- 
lations. 
3. Discussion of Results. 

The first approximation led to a characteristic equation that 
was linear ins and was very easy to inversely transform to find the 
time solution. The second approximation contained a characteristic 
equation that was quadratic in s and a little more difficult to trans- 
form, but still relatively simple. The quadratic term was a direct 
result of the 1/s factor in the second term of the series form of the 
Laplace Transform of v2. It follows that as succeeding terms are added 
to the approximation the characteristic equation will become of increasingly 
higher degree in s and more difficult to transform into the time domain. 
The requirement to solve difficult equations is contrary to the basic 
objective of finding a simple mathematical approximation. While the first 
approximation was exact in the steady-state and had the same general 
transient characteristics, the second approximation was very inaccurate, 


both in the steady-state and transient forms. If even more terms are 


16 


added in an attempt to find the best possible approximation, the 
resulting s domain equations become very difficult to solve and 


transform. 
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TIME 
.00 


ail 
.02 
“0S 
04 
AOS 
. 06 
cO/ 
.08 
Oe 


20 
soe 
40 
rod 


2.00 
S)5l018 
4.00 
elle, 
10.00 


Fore 


SOLUTION 


oo 


Ow 


EXACT 


. 00000 
.09828 
19656 
. 29483 
139310 
.49]37 
. 58963 
. 68788 
will 
.88435 
. 98256 
-JOs0u 
POOZ IO 
TOlooe 


4.88436 


24 


28. 


3] 


3G. 
£15 
36. 


.59594 
ay persile: 
77309 


8/7245 


.8]554 


05460 
40000 
40000 


APPROXIMATE 
SOLUTION 


> Ww PP 


eS 


AO 
24. 


ed 


38) 


34 
36 


- 00000 
.09815 
. 19603 
229300 
- 39100 
-48810 
. 98493 
. 68150 
.77781 
.8/7386 
.96965 
-91347 
. 74804 
./2636 


. 59674 


35320 
.56130 
63852 
36801 
.21503 
88677 
62811 
.40000 


TABLE III-1 


ERROR 


. 00000 
7900S 
.00053 
.00118 
.00210 
. 00327 
.00370 
. 00638 
.0083] 
.01049 
»0129) 
pUscZ 
.0841 1 
4421355) 
;2oroOg 
»24274 
2.19173 
3.73457 
4.50444 
4.60051 
Eo Maggers, 
V7 lies 
0.00000 











35} 


30 


ro) 


20| 


15 


10 


Oe 


time(sec. ) 


FIGURE III-] 
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B. APPLICATION TO A HOMOGENEOUS NON-LINEAR EQUATION 
The differential equation governing the discharge of a capacitor 
through a non-linear diode is given by Cunningham [ref. 3]. 


de 


C ape cle ar be = 0, e(0) = V 


= -Ae - Be, (111-6) 


where A = a/C and B = b/C. With the necessary condition that Gps Oz 
it can be seen from equation (III-6) that ep = 0. Recalling that the 


series form of the Laplace Transform of = is given by 


r @ @@ oo 


(en. 6) | (e. -aem) 
Lie} = £(s) | (eg ~ eg) - FO “EO. 


S 


Va) 


the terms in the series can be calculated. 


e+ Ae + Be“ = 0 Ep = -Ay - By< 
ee e e ee 66 = Z 2 2 3 e¢ e. 
e + Ae + 2Bee = 0 e. = AV + 3ABV™ + 2B°V~, er = 0 
Pees hee EES G0 ‘e) = -A°y-7ABV* - 10AB7V2 

- apry* - 2p°y4, e, = 0 


]. First Approximation. 


Ste*} = (e, 


Transforming equation (III-6) term by term, 


- €,)E(s) = -VE(s) 


Life + Ae + Bee = O} = sE(s) - V + AE(s) - BVE(s) = O (ITI-7) 
: V 
E(s) = s +A - BY 
Taking the inverse transform the approximate time solution is obtained 
e(t) = Vela eae (III-8) 


Cunningham [ref. 3] gives the true solution as 


Ve7At 
e(t) a Bir (III-9) 
1 - BV/A(e ~*-1) 


20 





Both solutions have initial values of V and final values of zero. 
The transient response is an exponential decay to 0 in both cases. 
A comparison between the two solutions is shown in Table III-2. It 
can be seen from the solutions that the approximation will be best 


for BV much less than A. 


2. Second Approximation. 


Including the first two terms of the series form in the 


approximation, the transform becomes 


(Semoue,,) 
_ (ep = @9) 
Sie“) = [lez - eg) - +] Els) 
After transforming equation (JII-6), 


; 2 
fe + Ae + Be“ = 0} = sE(s) - V + AE(s) + p (-y-SAVBY j= 0 


2,2 
E(s)(s + A - BV aC a = \ 


Vs 


oe aE 
s* + (A - BV)s - (ABV + B¢V*) 


Expanding in partial fractions, the s domain solution is found. 


M + 


ES) = 
5 + “BY , \a°-2By + aABy + Spey 


N 


~; (ABV) Re av + 4ABV SRV 


It is apparent that the denominator of the second term in the right-hand 
Side of the partial fraction expansion will, when transformed, cause the 
time solution to be unstable as the steady-state is approached. With the 
constants A= 1, V = 1, and B = .5, the time solution was found to be 


e(t) = .6052 exp(-1.15t) + .3946 exp(.75t) 


The second term is unstable. 


2] 


3. Discussion of Results. 

It can be shown that a third approximation to the solution 
using the first three terms in the series form of the transform of 
e*, with A= 1, B = .5, and V = 1, with result in the following s 
domain characteristic equation: 


2 


Eo) a — 
S) "5S: Bee 5S - 1),.5 


The corresponding time solution is 


.9965t , y.--45t 


e(t) = .22e K cos(1.2t + 6) 


As in the second approximation, this solution contains an unstable 
term. These results can be predicted by examining the characteristic 
equation in both cases. The presence of negative terms in these 
equations leads to unstable terms in the time solutions, since negative 
terms are contrary to known stability criteria. Furthermore, as more 
and more terms of the series approximation to the transform of the non- 


th order characteristic polynomial equation 


linear term are added, the n 
will contain at least one negative term if the (n-1) eh order equation 


contains at least one. 


C. DISCUSSION OF THE METHOD OF SUCCESSIVE APPROXIMATION 

In both examples of the use of Successive approximations, including 
only the first term in the approximation resulted in a satisfactory 
comparison to the exact solution while adding more terms not only made 
the solution more difficult, but also introduced unstable terms in one 
case, and an inaccurate solution in the other. In the following sections 
the final value method will be applied to more examples, but only the 


first approximation will be used. 


Ze 








EXACT APPROXIMATE 
TIME SOLUTION SOLUTION ERROR 
00 1.00000 1.00000 .00000 
Gi 98515 .99501 . 00986 
02 .97059 .99005 .01946 
.03 95631 98511 02880 
04 94232 .98020 .03782 
05 92859 97531 04672 
06 .91512 .97045 06467 
07 90191 96561 .06370 
08 88894 .96079 .07185 
.09 87622 .95600 .08918 
10 .86374 95123 .08749 
20 .75069 90484 .25415 
. 30 65583 .86071 .20488 
40 57546 .81873 24327 
50 50682 .77880 .27198 
1.00 .27953 60653 . 32700 
2.00 .09449 . 36788 27339 
3.00 .03375 22313 .18938 
4.00 .01229 .13534 12305 
5.00 .00450 08208 .07758 
10.00 .00003 .00708 .00705 
TABLE III-2 
= Tig B > Wicd Veet 


Zo 


IV. A DETAILED ANALYSIS OF THE APPLICATION OF THE FINAL 
VALUE METHOD USING FIRST APPROXIMATIONS ONLY. 
In the following sections the results obtained from the application 
of the approximation method to the equation 


Ye x + AX? 


= 1, x(0) = 0 (IV-1) 
will be shown. The effect of varying A, the coefficient of the non- 
linear term, will be studied by comparing the approximate solutions to 
solutions obtained on the SDS 9300 Digital Computer using Subroutine 
RKAM. In addition, the solutions obtained when a second derivative 


term was added to the same equation and A again varied will be analyzed 


in the same manner. 


A. APPLICATION TO A SECOND ORDER NON-LINEAR EQUATION 
To illustrate the application of the approximation method to equation 


(IV-1) consider this equation with A = 1: 
x +x + x° = 1, x(0) = 0 (IV-2) 


Transforming the equation term by tern, 


Six +x + x° = 1} = sX(s) - x(0) + X(s) + tx} = l/s 


1. Finding The Approximate Laplace Transform of ce, 


The approximate transform of the second order non-linear term, 
x, was found as follows: 
by definition, See = yes tut 
Expanding by integrating by parts, 


a 


[reat = xP | xe"Stat ~ (x2) [I xe S*dtdt + (x) | [xe SNatatat... 
0 oO ~% 
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Since othe = X(s), and integrating in time corresponds to 
multiplying by 1/s in the s domain, evaluating x2 and its derivatives 


at the extremes results in the series, 


Fix?) . (xé - x>)X(s) - (x2 - x2 XS) (xé - x ay we 


Rearranging terms, 


Co) 
Lix®} = X(s) | 02 - xr) - i ee - | (1V-5) 


For a first approximation all derivatives of x2 are neglected and the 


linearized approximation to the Laplace Transform of 3 becomes 


Sixth = (x6 = x4)X(s) (1V-6) 


Recalling the required condition that Xe = 0, the final value, Xe» can 


be found. 
3 


{f 
—! 


Cat ar 


3 +x - | 


i} 
© 


(IV-7) 


The cubic equation has three roots, one real and two complex. Selecting 
the real one, 
Xe = .682 
Using this value of Xe» equation (IV-6) becomes 
Lix?} = .4651X(s). 
2. Alternate Methods of Finding The Approximate Transform. 

To further illustrate the obtaining of the approximate transform 
of the non-linear term, two different approaches will be used. Equation 
(IV-3) can also be expanded as 

(£0 


oR >, | : 
| woe Scat = xf xPerstet - | | ule Stutdy nan | Nae Stutdtdtee+: 
o i Ly 


ee, 


Z3 


Using the Laplace Transform of x2 found by Brady [ref. 2], 
x7e Stat = sx@(s), 


and evaluating the terms in the expansion as before, 


® e ee 


Fix} = (X¢ - x9 )8X“(s) - ~~ Xo! sX*(s) + OF a a 0) sX2(g)-0 5: 
Neglecting derivatives and rearranging, 
Lix?} = sX"(s)x, (IV-8) 
Now, expanding X(s) in its infinite series form, 
X(s) = SE cullen (ep - %o) 4 Op = Xo) is 
S 2 3 


S S 
If derivatives are again neglected, 


Xe - X x 
ee 
X(s) = S S 
and equation (IV-8) becomes 
SX -x-X(S) 
an - tt = x2X(s) 


The approximate transform thus obtained is identical to equation (IV-6). 
The final approach will be to take the transform of x? directly, 


using results found by Brady [ref. 2]: 


ix} = 5°x3 (5) = 5°X(s)X“(s) ' (IV-9) 


Writing x*(s) as the product of two infinite series, 


mpieialmtipctohpl Opr%) 70) i jeiarae’s ‘aa - 


S 


Neglecting derivatives and products containing derivatives, 
2 

ee ee 

X"(s) = S ser 














Substituting this result into equation (IV-9), the approximation becomes, 


Z 2 
s“X(s)x 
Six?) = 5 = xEX(s) 
S 


Again, the approximate transform obtained is identical to equation 
(IV-6). 
3. Using The Approximate Transform To Find The Solution. 
Completing the term by term transform of equation (IV-1) 
: Ss a 2 " 
Six +x+x> = 1} = sx(s) - x(0) + X(s) + X'eX(s) = ]/s 
With x(0) = 0, this equation becomes 


SAS NAG alls 


_ ] 
X(s) = S(s + 1.4651) 


Taking the inverse transform, the approximate solution is found. 


[> 


: A. B 
A(s) ““Seeame) aa. s ~ (eee CELE 


. ] ‘ -] 
e-seao >  ”— aaledesin 


_ 6825 6825 
Bee eT AGS 


-1.4651t 


x(t) = .6825 - .6825e 


A comparison between the approximate solution and the computer solution 
is shown in Table IV-3. 
4. Variation of the Coefficient of the Non-linear Term. 
Extending the methods developed to the solution of the general 
equation, 
eal lt Ax® Sallie OD SeOx 


the final values are found by solving the cubic equation, 


ai 


With A = .5, this equation becomes 


Coy 2 ener 


The solution to the cubic equation of the form, 


x2 + ax + b= 0 


can be found [ref. 4]: 


Letting A = \- of ee, B = \-2 - - So, 


then the roots are 


{ae Ly 2 et, : ia) - AB) ES 


Choosing the real root, x =A+B. 


With A = .5, the auxiliary equations are: 


_\i2 a44 , 8 owe. oe =e 
Me (eer , e+ \3 n+ 


= —_ 
¢ = +7713, x, = 5949. 


Similarly, for other values of A, the values of Xe and xe were found 


and x 


to be as follows: 


2 


A *f ae 
5 STS ~5949 
9 .6972 4861 
jet, ~6825 .465] 
1.1 .6709 ~4500 
peas ~6257 3914 
2.0 .5904 348] 
owed 0, 4847 2349 


The general characteristic equation was found by transforming the general 


equation term by term, with the initial condition, x(0) = 0. 
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sae + x + Ax? = ]} = sX(s) + X{s) + Ax6X(s) = l/s 


| oO 


C 
X(s) — = Se 
Xie 1+ Axe) (sg + ] + Axe) 


] c -| 


oo 9 —= 9 
1 + Axe 1+ Axe 


The generalized time solution (approximate) is then 


Z 
mt) = 2. a y) aU. uy AX'g)t 
He: AX¢ Ie AXe 





With the same values of A, the time solutions were found to be as 





follows: 

A x(t) 

5 .7708 - .7708 exp (-1.2975t) 
9 6905 - .6956 exp (-1.4375t) 
me 6825 - .6825 exp (-1.4651t) 
All .6689 - .6689 exp (-1.4950t) 
ilies .6300 - .6300 exp (-1.5871t) 
2.0 5896 - .5896 exp (-1.6962t) 
5.0 .4600 - .4600 exp (-2.1745t) 


Tables IV-] through IV-7 show comparisons between the approximate 
solutions and the solutions found on the digital computer for the 
several values of A. For all values of A the approximations were 

quite satisfactory. In each case the approximate solution rose slightly 
Slower to the steady-state. For the worst case, A = .5, the approximate 
Steady-state value differed from the computer solution bt 2.6%. The 
maximum difference in the transient values was approximately 10%, but 


typical differences were much less. 
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TIME 


KSEC) 


oS] Fo Toso 2 2&2 Go G2 © 


Oo OO © 


a, S&S = Se oweoe GF Oo OD 


| De 


Oo Oo SP W 


~l 


COMPUTER 
SOLUTION 


. 00000 
.09515 
le) 
120801 
.32784 
. 38958 
.44418 
-49215 
.5340] 
- 57030 
.60158 
. 74303 
./6665 
102) 
;77U82 
./7090 
re 7 OOM 
eT TUIL 
OI 
sive 


APPROXIMATE 
SOLUTION 


.00000 
-02570 
- 17618 
. 24853 
.31209 
- 36790 


-41693 


-45999 


-4978] 
oe Os 
.96021 
wa oOZ2 


DUG 
./6651 
O50 
./7048 
-/7071 
-77078 
.//079 
rT OLG 


TABLE IV-] 


A=0.5 


30 


ERROR 

00000 
.00136 
.00494 
00998 
.01575 
.02168 
02725 
.03216 
.03620 
03927 
.04137 
02941 
.01157 
00376 
.00118 
00042 
00021 
.00014 
.00013 
.00012 





TIME 


(SEC) 


© 
LO 


Oo © 


COMPUTER 
SOLUTION 


.00000 
.09514 
. 18100 
-25/98 
. 32640 
. 38656 
43884 
-48376 
,52193 
.55402 
socw75 
.68374 
.69499 
.69612 
168623 
.69624 
.69624 
.69624 
.69624 


APPROXIMATE 
SOLUTION 


TABLE IV- 


A= 0.9 


3] 


.00000 
.09314 
.17381] 
. 24367 
. 30418 
SIDS, 
.40198 
4430 
»47535 
50484 
. 53038 
.65636 
.68628 
BOOS 
-69507 
369538 
oo oe, 
.69559 
soo 560 


Z 


ERROR 


. 00000 
.00200 
.00719 
.0143] 
sO2i22 
NOR: 
- 03686 
.04246 
.04658 
.04918 
.05037 
.02738 
.0087 1 
200273 
.00116 
.00076 
. 00067 
.00065 
.00064 


TIME 
(SEC) 


2) fen Sk Se 8 Se a <> Re <= ee > Re an 
Oo oF: FF WW 


a / 
Co O™ 


Ho O98 FSP WD PO 
Se eas S-lCOCO]lClC COCO 


COMPUTER 
SOLUTION 


. 00000 
.09514 
mie097 


#25785 
.32605 


#88583 
-43756 
.48176 
ATTT 
seoUzZ2 
no7076 
.67181 
.68138 
.68224 
.68233 
HOG233 
.68233 
.68233 


APPROXIMATE 
SOLUTION 


TABLE IV 
A = 1.0 


SZ 


. 00000 
-09301 
a 335 
.24274 


. 30267 
35443 
.39914 


-43776 


47111 
-49992 


-5248) 
. 64606 
.67408 
. 68005 
.68205 
. 68240 
.68248 
.68249 


~3 


ERROR 


.00000 
.00213 
.00762 
.01511] 
mz 500 
.03140 
.03842 
.04400 
04796 
.05030 
8051/15 
102575 
.00730 
.00169 
.00028 
-00007 
.00015 
.00016 








TIME COMPUTER APPROXIMATE 

(SEC) SOLUTION SOLUTION ERROR 
0.0 00000 00000 .00000 
0.1 09513 09288 02259 
0.2 .18094 17287 00807 
0.3 25772 ,24175 01597 
0.4 32569 30106 02453 
0.5 38508 35214 03294 
0.6 43626 39613 .04013 
0.7 47976 43401 04575 
0.8 51626 46662 04964 
0.9 54650 .4947) ,05179 
1.0 57131 51890 05241 
2.0 66073 63526 02547 
3.0 66891 66136 00755 
4.0 66960 66721 00239 
5.0 66966 66852 00114 
6.0 66966 66881 00085 
7.0 66966 66888 .00078 
8.0 66966 66890 00076 

TABLE IV-4 
A= 1.) 
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APPROXIMATE 


TIME COMPUTER 

(SEC) SOLUTION SOLUTION ERROR 
0.0 .00000 .00000 00000 
0.1 .09513 09246 00267 
ee . 18082 .17134 00948 
0.3 , 25720 . 23865 .01855 
0.4 32438 . 29609 .02819 
0.5 38218 34509 .03709 
0.6 43126 .38690 04436 
0.7 47211 42258 .04953 
0.8 .50556 .45302 05254 
0.9 53256 .47899 .05356 
1.0 .55406 .50115 .05291 
2.0 62321 .60365 .01956 
3.0 .62787 62461 00326 
4.0 62816 .62890 .00070 
5.0 .62818 62877 .00159 
6.0 62818 62995 .00177 
7.0 62818 62999 .00181 
8. 62818 .63000 .00182 

TABLE IV-5 
A= 1.5 
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TIME COMPUTER APPROXIMATE 

(SEC) SOLUTION SOLUTION ERROR 
0.0 00000 .00000 00000 
0.1 .09511 .09199 .03129 
0.2 .18067 16962 .01105 
0.3 25655 .23514 02141 
0. 32255 .29044 02311 
0.5 37867 .33712 04155 
0.6 42529 .37651 04878 
0.7 .46315 .40875 05340 
0.8 .49329 .43781 05548 
0.9 .51687 .46149 05538 
1.0 53506 .48148 04356 
2.0 .58710 56977 .01733 
30 58963 58596 00367 
4.0 58975 58893 00082 
5.0 .58975 .58948 00027 
6.0 58975 58958 00027 

TABLE IV-6 
A= 200 
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TIME 


(SEC) 


Oo OO & 


TOR 


Oo OF © 


oT > FSF W 


oT Oo oO 


COMPUTER 
SOLUTION 


00000 
.09505 
17977 
25276 
.31280 
35978 
39484 
.41999 
43748 
44938 
45734 
.47075 
4723) 
47249 
4725) 
47251 
4725) 
4725) 
4725) 
4725) 
47251 


APPROXIMATE 
SOLUTION 


TABLE IV- 


A = 5.0 


36 


- 00000 
.08990 
Bi6223 
~22042 
.26724 
. 3049] 
Bo 3022 
#9595] 
ao o23 
2250] 
-40771 
44237 
.45406 
.45800 
- 45932 
-45977 
.45992 
45997 
ogg 
- 46000 
-46000 


/ 


ERROR 


. 00000 
.00515 
.01754 
. 03234 
.04556 
.05487 
lloySley4 
. 06038 
"Wo0c5 
.05437 
.04963 
.02838 
#00825 
.01449 
.01319 
.01274 
OT 59 
.01254 
.01252 
<025] 
eo25 | 


B. THE EFFECT OF ADDING A SECOND DERIVATIVE TERM TO THE NON-LINEAR 
EQUATION 


A second derivative term was added to equation (IV-1) and the 
approximate transform was used to find an approximate solution to the 
resulting equation. The equation, 


ee Ae = ei) 0, (IV-11) 


still contained a second-order non-linear term, x, 


1. Use of the Final Value Method to Find the Approximate Solution. 
As in Part A the solution method will be illustrated on the 


example equation with A = 1: 





a  eeOluae (I-12) 


The final value is the same as the one found for equation (IV-1), since 
Xe will also be zero in the steady-state. The approximate transform 
of the second-order non-linear term is also the same since its derivation 
was independent of the system equation. 

Transforming equation (IV-12) term by term, which x(0) = 0, 


3 2 


= ]} = 5°X(s) + sX(s) + X(s) + Xe 


S Ne xe x X(s) = l/s 


The § domain characteristic equation is then 


] 
+ s + 1.4651) 


X(s) = 


5(5° 


Expanding in partial fractions the s domain solution can be found. 





oe ele ae A B ‘i Bx 
s(s° +s +1.4651) °* ° + 5 - j1.102 © 5 + .5 + j1.102 

A = .6825 

o 
ae 23! 
B = -.3600e 92" 
=J24 23" nicdee 3) 
ig) = ESE SoS -.36¢e 


5 50 8] Se 2 Grae mI 


oy 





The corresponding time solution was found to be 


fois 


x(t) = .6825 - .7200e— cos (1.102t - .42547) (IV-13) 


Equation (IV-13) is basically that of a second order system, which was 
expected since the example equation is basically second order, with a 
non-linear term. A satisfactory approximation should start at zero and 
follow a damped sinusoidal path to some steady-state value. The approxi- 
mation here does have a damped sinusoidal characteristic and proceeds to 
a steady-state value that is the same as the value predicted. To meet 
the requirement that the solution begin at zero it is necessary that 


the second term in the approximate solution be equal to zero, or that 


.72 cos (-.42547) = .6825 

It will be shown that this is very nearly the case. 

2. Solution To The General Equation. 

The characteristic equation for the general case, equation 

(IV-11), was found to be 

] 
ee ee Axé) 
For the same values of A as in Part A the time solutions were found to 


be as follows: 





5 7708 - .8468e “2tcos (1.024t - .45378) 
9 6956 - .7642e “Pecos (1.091t - .42980) 
.0 6825 - .7200e *2rcos (1.102t - .42547) 
] .6689 - .7526e “5 ECOS (1.116t - .41384) 
5 .6300 - .6780e “prcos (1.155t - .40833) 
0 .5896 - .6400e “’cos (1.201t - .42927) 
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For all values of A the approximate solution oscillated slightly 
Slower than the computer solution, but both solutions had the same 
number of oscillations. The difference between peak overshoots was 
in all cases much less than 10%, as was the difference between final 
values. As can be Seen in Tables IV-8 through IV-13, the initial 
values of the approximate solutions were indeed very close to zero, 
the worst case being when A = 1.0 and A = 1.1, when the value of 


the approximate solution at t = 0 was -.00217. 


Shs, 


TIME COMPUTER APPROXIMATE 
SEO SOLUTION SOLUTION 
0.0 .00000 .00097 
= 10440 11243 
oe : 33989 » 33845 
1.5 .6045] -58305 
24:0 iB lee .77810 
ane) . 92862 .89458 
2.8 .94273* . 92666 
ial .92556 .93484* 
BaD 5O7 30] 91794 
4.0 .19738 .87134 
4.5 . 74154 .81807 
5.0 .71795 .7740) 
5.2 .71687# . 76086 
6.0 .74152 .13543# 
65 . 76366 . 73807 
J 20 . 77966 .74757 
A: .78646* .76130 
O70 . 78467 ./690] 
8.5  /576 .77546 
9.2 .77029 .77844* 
18 ./6646 ./7660 
TABLE IV-8 
A=0.5 
* maximum 
# minimum 
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ERROR 


.00097 
.00803 
.00144 
.02146 
03946 
.03404 
.02607 
00928 
.04463 
.07396 
.07653 
.05706 
.04399 
00509 
02559 
.03209 
.02516 
01566 
.00332 
.00815 
.01014 








TIME COMPUTER APPROXIMATE 

(SEC) SOLUTION SOLUTION ERROR 
0.0 .00000 .00090 .00090 
0.5 . 10440 . 10442 .00002 
2 . 33956 - 32977 .00979 
| .59980 S675 .03275 
2.0 . 7987 2 . 74632 .04740 
2.5 .86480* .84110 .02370 
Ps. 83921 .86018* .02097 
$.5 .7395] .82436 .08485 
4.0 .66772 . /6820 .09048 
4.6 .63470# 505 .07035 
5.0 .64183 .67630 03447 
5.2 .68537 .65663# .02874 
6.5 70275 . 66807 .04468 
6.9 .71623* .67896 .03727 
sees .70913 .69378 . 09535 
8.0 .69914 . 70162 .00221 
8.6 .69123 .70484* .01361 
9.0 .68950# . 70408 .01458 
9.9 69827 69876 .00549 

TABLE IV-9 
A = 0.9 
* maximum 
# minimum 


4) 


TIME COMPUTER APPROXIMATE 
SEC) SOLUTION SOLUTION 
0.0 .00000 -.00217 
0.5 . 10439 . 10442 
1.0 . 33947 .32977 
0 . 59862 . 56705 
2.0 78797 . 74632 
255 .85002* .84110 
2.9 .81757 .86018* 
255 A SL5 .82436 
4.0 .64515 . 76820 
4.5 .61960# .71417 
5.0 .63091 .67630 
are .67750 .65663# 
6.7 10314 # .67342 
inte .69274 .69378 
8.0 "oe2o7 .70162 
8.6 .67578 . 70484* 
9.0 67525 . 70408 
9.9 .68057 .69876 

TABLE IV-10 

A = 1.0 

* maximum 
# minimum 
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ERROR 
00217 
00003 
.00970 
03157 
.04165 
.00912 
04261 
.11061 
12305 
09457 
04539 
02087 
02972 
.00104 
01895 
02906 
.02883 
.01819 





TIME COMPUTER APPROXIMATE 
(SEC) SOLUTION SOLUTION ERROR 
0.0 .00000 -.00271 .00271 
@.5 . 10440 .08885 O55 
1.0 . 33939 . 32040 .01899 
as. .59747 .56024 03723 
ZnO .7823] . 73669 .04562 
z.5 .83615* 82438 OM 7 
ass .83615 82438 .01177 
2.8 .81156 .Gavo4* .02598 
B.5 .69032 os . 10141 
4.0 .62529 oS .11324 
4.5 .60557# OFIMS .07158 
5 al) .62180 .64182 .02002 

.6 .65910 .62765# .03146 
6.0 .6794] .63143 .04798 
6.6 69RzZ7* .64712 .04415 
ee .68784 . 65890 .02894 
Teds .67746 .67071 .00674 
8.0 .66750 <67 735 .00015 
8.6 66214 .67882* .01668 
a. / .66202# .67853 .01831 
9.5 .66608 .67360 §0C7 25 
9.9 .66911 .67078# .00167 

TABLE IV-11 
A = 1.1 

* maximum 
# minimum 
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TIME COMPUTER APPROXIMATE 
(Sa SOLUTION SOLUTION 
0.0 .00000 .00019 
G5 . 10440 - 10951] 
0 . 33906 . 32818 
5 . 59286 .5518] 
2.0 . 76043 er inO3 
253 SIS . 76455 
Zo . 76030 . 78966* 
365 .61477 m9 | 
GisZ .56072# .65212 
a0 .59678 .59613 
5.5 .63090 .58906# 
OZ soe! * .60276 
G20 .63870 .62639 
Tes FO2Z622 .6362] 
ear .61907# .64051* 
9.0 .62442 oie / 
S08 .62902 .63245 
9.9 .63112 .62985# 
TABLE IV-12 
* maximum 


SH: 


minimum 


ERROR 


.00019 
.00511 
.01088 
.04105 
.04940 
.02696 
.02936 
- 11904 
00860 
.00065 
.04184 
.04945 
e023 | 
. 00999 
.02144 
gil) sis 
.00343 
.00127 





a e er rere 


TIME 
(SEC) 


on 


oro @& Bei ®t 


He 


“a NM OD 


OO © 
© 


maximum 
minimum 


COMPUTER 
SOLUTION 


.00000 
.10439 
33865 
58720 
.74929* 
.71288 
63342 


.54876 


519454 
53894 
57814 
59863 
.59863* 
60512 
579214 
57947 
58527 
59100 
.59385* 


APPROXIMATE 
SOLUTION 


-00011 
.09845 
31 138 
. 52996 


11719 
.74679* 


.71233 
67929 
61830 
57207 
54967 
547084 
5548) 
57821 
60105 
60111 
59816 
59371 
58848# 


TABLE IV-13 


ERROR 


.00011] 
.00594 
.02727 
.05724 
.03210 
.0339] 
.0789] 
. 13058 
. 09885 
 OSSWS 
.07047 
.05155 
.05769 
.02691 
.02184 
.02611] 
.01289 
.0027] 
00537 


V. CONCLUSION 


The specific object of the methods developed was to find a 
mathematically simple procedure with which to approximate the 
solution to certain non-linear differential equations. The only 
requirement was that the steady-state (final) value of the system 
variable could be obtained. This requirement implies a non- 
oscillatory steady-state, where all derivatives are zero. The 
best results were obtained in the case where there was an applied 
forcing function (non-homogeneous equation). The nature of the 
solution does not limit this function to a unit step, although 
either a constant input or one which decays to zero in the steady- 
State are the only inputs which will permit a final value of the 
system variable to be found. For the discharging capacitor example, 
with no input applied, the transient response was not very accurate, 
although the general path followed and the final value obtained were 
essentially the same as for the exact solution. 

The usefulness of the final value method of approximation lies 
in its ready use by the engineer. The mathematics required is minimal 
and the solutions obtainable are sufficiently close to the exact solution 
to allow its use in analyzing non-linear systems which have a constant 


Steady-state. 
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TIME COMPUTER APPROXIMATE 
(SEC) SOLUTION SOLUTION ERROR 
0.0 . 00000 .0001 1 .0001 1 
.5 . 10439 .09845 .00594 
1.0 . 33005 Loe .02727 
Ps Dove” s52e06 .05724 
Lop .74929* .71719 05270 
2.6 .71288 .74679* .03391 
Stoll .63342 m2 38 .0789] 
325 54876 .67929 . 13053 
4.0 519454 .61830 09885 
4.5 53894 107207 .03313 
5.0 .57814 . 54967 .07047 
5.3 59863 .54708# 05155 
5.9 J59663* 5548] 105769 
6.5 S605 1/2 57871 .02691 
7.8 579214 60105 02184 
: 7.9 57947 .60111* 02611 
8.5 -5652/ .59816 .01289 
9.0 .59100 59371 .00271 
9.7 .59385* 588484 .00537 
TABLE IV-13 
A = 2.0 
* maximum 
# minimum 





V. CONCLUSION 


The specific object of the methods developed was to find a 
mathematically simple procedure with which to approximate the 
solution to certain non-linear differential equations. The only 
requirement was that the steady-state (final) value of the system 
variable could be obtained. This requirement implies a non- 
oscillatory steady-state, where all derivatives are zero. The 
best results were obtained in the case where there was an applied 
forcing function (non-homogeneous equation). The nature of the 
solution does not limit this function to a unit step, although 
either a constant input or one which decays to zero in the steady- 
state are the only inputs which will permit a final value of the 
system variable to be found. For the discharging capacitor example, 
with no input applied, the transient response was not very accurate, 
although the general path followed and the final value obtained were 
essentially the same as for the exact solution. 

The usefulness of the final value method of approximation lies 
in its ready use by the engineer. The mathematics required is minimal 
and the solutions obtainable are sufficiently close to the exact solution 
to allow its use in analyzing non-linear systems which have a constant 


steady-state. 
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